Molecular Docking Analysis of Triazole Analogues as Inhibitors of Human Neutrophil Elastase (HNE), Matrix Metalloproteinase (MMP 2 and MMP 9) and Tyrosinase

 

V. Vijayakumar1, N. Radhakrishnan2*, P. Vasantha-Srinivasan3

1Department of Chemistry, Anna University, Chennai-25, Tamil Nadu, India.

2,3Department of Biochemistry, St. Peter’s Institute Higher Education and Research, Avadi, Chennai,

Tamil Nadu, India.

*Corresponding Author E-mail: nrkishnan@gmail.com

 

ABSTRACT:

Triazole analogues are well known to have desirable properties for medicinal chemistry such as hydrogen bonding capability, moderate dipole character, rigidity and stability under in vivo conditions etc. This prompts us to carry out the present study, we have selected 20 triazole analogues and evaluated on the docking behavior of four targeted enzymes such as HNE, MMP 2, MMP 9 and tyrosinase. The molecular physicochemical analysis revealed that all the tested ligands (expect ligand 5 and 18) showed nil violation and complied well with the Lipinski’s rule of five. ADME analysis showed that ligand 15 alone predicated to have low gastrointestinal (GI) absorption effect. Docking studies revealed that ligand 18 [2-(2,4-difluorophenyl)-1-(4-((1-(4-ethylbenzyl)-1H-1,2,3-triazol-4-yl)methyl)piperazin-1-yl)-3-(1H-1,2,4-triazol-1-yl)propan-2-ol] had exhibited the maximum binding energy for three target enzymes such as HNE, MMP 2 and MMP 9. The present study has paved a new insight in understanding 20 triazole analogues as potential inhibitors against HNE, MMP 2, MMP 9 and tyrosinase.

 

KEYWORDS: Molecular Docking, Triazole derivatives, Human neutrophil elastase, Matrix metalloproteinase and Tyrosinase.

 

 


INTRODUCTION:

Triazole analogues are five member heterocyclic compounds basically having two carbon and three nitrogen atoms in their structure. Triazole analogues have been reported to posses various biological activities such as analgesic, analeptic, anti-anxiety, anti-cancer, anti-convulsant, anti-depressant, anti-fungal, anti-histaminic, anti-inflammatory, anti-hypertensive, anti-malarial, anti-microbial, anti-oxidants, anti-tubercular and anti-viral etc1.

 

Moreover, some commercial available drugs like Alprazolam (anti-anxiety drug), Anastrozole (breast cancer), Etizolam (anti-anxiety), Fluconazole (anti-fungal), Intraconazole (anti-fungal), Letrozole (breast cancer), Ribavirin (anti-viral), Rizatriptan (anti-migraine), Rufinamide (anti-convulsant), Triazolam and Voriconazole (anti-fungal) are known to have triazole moiety2. Apart from these, triazole analogues have been reported to posses unique advantages like i) potent pharmacological activities, ii) no or low adverse effects, iii) more bioavailability, iv) poor drug resistances and v) better pharmacokinetics property etc.1.

 

Molecular docking is one of the bioinformatics tools to study the interaction between the small molecules (ligands) and protein/enzyme (receptor), predict the binding site of the ligand. Previously, researchers had reported the molecular docking studies of (2E)-N-Hydroxy-3-[3-(Phenylsulfamoyl) Phenyl] prop-2-Enamide (Belinostat) against Ebola virus glycoprotein and 1, 2, 4-triazin analogue of diclofenac as potential ligand for Parkinson’s disease2-3. Panigrah and co-workers had reported the synthesis, antimicrobial activity and molecular docking studies of novel oxazolidinone-thiophene chalcone hybried derivatives4. Chelamalla and Makula reported the molecular docking and ADMET analysis of pyrimidine coumarin scaffolds5. Hasan and coworkers had reported the molecular docking investigations of 3, 5, 7-trihydroxy-6,7-dimethoxy flavone6. Ramjith and Shahim had reported the molecular docking studies of novel imidazo[2,1-b]-1,3,4 thiadiazole derivatives7.  Similarly, few other organic compounds have been reported to possess biological activity both in vitro and in vivo studies8-10.

 

Kumar and co-workers has reported anti-inflammatory activity for 28 (1, 2, 4)-triazole derivatives, among which triazole derivative 17 (n-butyl group at position 4) had shown better in vivo anti-inflammatory activity11. Similarly, Wuest and   co-workers has reported that diaryl substituted 1, 2, 3 triazole derivatives have exhibited good cyclooxygenase (COX-2) inhibition activity12. This prompts us to carry out the present study, we have selected 20 triazole analogues they are (1) 4H-1,2,4-triazole; (2) 4-butyl-4H-1,2,4-triazole; (3) 4-octyl-4H-1,2,4-triazole; (4) 4-(4H-1,2,4-triazol-4-yl)butan-1-amine; (5) 4-octyl-3,5-dipropyl-4H-1,2,4-triazole; (6) 4-amino-5-ethyl-4H-1,2,4-triazole-3-thiol; (7) 3-ethyl-5,6-dihydro-7H-[1,2,4]triazolo[3,4-b][1,3,4]thiadiazin-7-one; (8) 5-methyl-4-phenyl-4H-1,2,4-triazole-3-thiol; (9) (E)-2-(((3-(p-tolyloxy)-4H-1,2,4-triazol-4-yl)imino)methyl)phenol; (10) 2-(4-acetylphenyl)-5-methyl-2H-1,2,3-triazole-4-carboxylic acid; (11) 5-azido-4H-1,2,4-triazol-3-amine; (12) 2-(2,4-difluorophenyl)-1,3-di(4H-1,2,4-triazol-4-yl)propan-2-ol; (13) 4-(2-(2,4-dimethylphenyl)-2-methyl-3-(1-nitro-1,5-dihydro-4H-1,2,4-triazol-4-yl)propyl)-4H-1,2,4-triazole; (14) 6-((R)-3-aminopiperidin-1-yl)-3-((1-methyl-4,5-dihydro-1H-1,2,3-triazol-4-yl)methyl)-1-(prop-1-yn-1-yl)dihydropyrimidine-2,4(1H,3H)-dione; (15) (2-(4-((2,4-dioxo-3,4-dihydropyrimidin-1(2H)-yl)methyl)-2H-1,2,3-triazol-2-yl)ethyl)phosphonic acid; (16) (2R,3S,4R,5R)-2-(hydroxymethyl)-5-(5-(trifluoromethyl)-1H-1,2,3-triazol-1-yl)tetrahydrofuran-3,4-diol; (17) 1-(1-(2,4-dichlorophenyl)-2-(1H-1,2,3-triazol-1-yl)ethyl)-4-ethyl-5-methyl-1H-1,2,3-triazole; (18) 2-(2,4-difluorophenyl)-1-(4-((1-(4-ethylbenzyl)-1H-1,2,3-triazol-4-yl)methyl)piperazin-1-yl)-3-(1H-1,2,4-triazol-1-yl)propan-2-ol; (19) 4-((2R)-3-(2,4-difluorophenyl)-4-(1H-1,2,4-triazol-1-yl)butan-2-yl)pyrimidine; (20) 4-(2-((2S,3R)-3-(2,4-difluorophenyl)-3-hydroxy-4-(1H-1,2,4-triazol-1-yl)butan-2-yl)-2,3-dihydrothiazol-4-yl)benzonitrile, were evaluated on the docking behavior of Human neutrophil elastase (HNE), matrix metalloproteinase (MMP 2 and MMP 9) and tyrosinase. The results shows that the potential of triazolium analogues as human neutrophil elastase (HNE), matrix metalloproteinase (MMP 2 and MMP 9) and tyrosinase inhibitors.

 

MATERIALS AND METHODS:

Ligand preparation:

Chemical structures of 20 ligands (triazole analogues) were drawn in Chem Bio Draw Ultra 12.0 and energy minimization of ligands was carried out by Chem Bio 3D Ultra 12.0, according to the reported procedure13. These structures were employed for further Patch Dock.

 

Target protein identification and preparation:

The three dimensional (3D) structures of the HNE (PDB ID: 1H1B with resolution of 2.00 A) and MMP 2 (PDB ID: 1QIB with resolution of 2.80 A), MMP 9 (PDB ID: 4H1Q with resolution of 1.59 A) and Tyrosinase (PDB ID: 2Y9W with resolution of 2.30 A) were obtained from the Research Collaborator for Structural Bioinformatics (RCSB) Protein Data Bank (www.rcsb.org)13-14. A chain of these proteins was processed individually by removing another chain (B, C and D), ligands in addition to the crystallographically observed water particles (water without hydrogen bonds). The protein mentioned above was prepared using UCSF Chimera software (www.cgi.ucsf.edu/chimera).

 

ADME analysis:

ADME (Absorption, Distribution, Metabolism and Excretion) analysis was performed by Swiss ADME analysis was carried by a standard default protocol15.

 

Docking studies:

Docking studies were carried out by the PatchDock online server (http://bioinfo3d.cs.tau.ac.il/PatchDock). PatchDock adopts geometry-based molecular docking algorithm method was used to recognize the binding scores, by binding residues atomic contact energy of the given ligands. The docking results were obtained through the email address. We also use to get uniform resource locator (URL) which provides the top 20 solutions in a table form via email. From these, the top one solutions (the docked protein-ligand complex) was selected and downloaded in a database (pdb) file format. Further, the binding site analysis was carried out by PyMOL software (www.pymol.org)16.

 

RESULTS AND DISCUSSIONS:

Twenty triazole analogues have been selected for the present study, as shown in the Table 1.

 


Table 1: Represents the International Union of Pure and Applied Chemistry (IUPAC) name and Simplified Molecular Input Line Entry System (SMILES) of 20 triazole analogues.

S. No

IUPAC name

SMILES

1

4H-1,2,4-triazole

C1=NN=CN1

2

4-butyl-4H-1,2,4-triazole

CCCCN1C=NN=C1

3

4-octyl-4H-1,2,4-triazole

CCCCCCCCN1C=NN=C1

4

4-(4H-1,2,4-triazol-4-yl)butan-1-amine

NCCCCN1C=NN=C1

5

4-octyl-3,5-dipropyl-4H-1,2,4-triazole

CCCC1=NN=C(CCC)N1CCCCCCCC

6

4-amino-5-ethyl-4H-1,2,4-triazole-3-thiol

NN1C(CC)=NN=C1S

7

3-ethyl-5,6-dihydro-7H-[1,2,4]triazolo[3,4-b][1,3,4]thiadiazin-7-one

CCC1=NN=C2N1NCC(S2)=O

8

5-methyl-4-phenyl-4H-1,2,4-triazole-3-thiol

CC1=NN=C(S)N1C2=CC=CC=C2

9

(E)-2-(((3-(p-tolyloxy)-4H-1,2,4-triazol-4-yl)imino)methyl)phenol

OC1=CC=CC=C1/C=N/N2C(OC3=CC=C(C)C=C3)=NN=C2

10

2-(4-acetylphenyl)-5-methyl-2H-1,2,3-triazole-4-carboxylic acid

CC1=NN(C2=CC=C(C(C)=O)C=C2)N=C1C(O)=O

11

5-azido-4H-1,2,4-triazol-3-amine

NC1=NN=C(N=[N+]=[N-])N1

12

2-(2,4-difluorophenyl)-1,3-di(4H-1,2,4-triazol-4-yl)propan-2-ol

OC(CN1C=NN=C1)(C2=C(F)C=C(F)C=C2)CN3C=NN=C3

13

4-(2-(2,4-dimethylphenyl)-2-methyl-3-(1-nitro-1,5-dihydro-4H-1,2,4-triazol-4-yl)propyl)-4H-1,2,4-triazole

CC(CN1C=NN=C1)(C2=CC=C(C)C=C2C)CN3C=NN([N+]([O-])=O)C3

14

6-((R)-3-aminopiperidin-1-yl)-3-((1-methyl-4,5-dihydro-1H-1,2,3-triazol-4-yl)methyl)-1-(prop-1-yn-1-yl)dihydropyrimidine-2,4(1H,3H)-dione

CN1CC(CN2C(CC(N3CCC[C@@H](N)C3)N(C#CC)C2=O)=O)N=N1

15

(2-(4-((2,4-dioxo-3,4-dihydropyrimidin-1(2H)-yl)methyl)-2H-1,2,3-triazol-2-yl)ethyl)phosphonic acid

O=P(O)(O)CCN1N=CC(CN(C(N2)=O)C=CC2=O)=N1

16

(2R,3S,4R,5R)-2-(hydroxymethyl)-5-(5-(trifluoromethyl)-1H-1,2,3-triazol-1-yl)tetrahydrofuran-3,4-diol

OC[C@@H]1[C@@H](O)[C@@H](O)[C@H](N2C(C(F)(F)F)=CN=N2)O1

17

1-(1-(2,4-dichlorophenyl)-2-(1H-1,2,3-triazol-1-yl)ethyl)-4-ethyl-5-methyl-1H-1,2,3-triazole

ClC1=CC=C(C(N2N=NC(CC)=C2C)CN3C=CN=N3)C(Cl)=C1

18

2-(2,4-difluorophenyl)-1-(4-((1-(4-ethylbenzyl)-1H-1,2,3-triazol-4-yl)methyl)piperazin-1-yl)-3-(1H-1,2,4-triazol-1-yl)propan-2-ol

FC1=CC(F)=C(C(CN2C=NC=N2)(O)CN3CCN(CC4=CN(CC5=CC=C(CC)C=C5)N=N4)CC3)C=C1

19

4-((2R)-3-(2,4-difluorophenyl)-4-(1H-1,2,4-triazol-1-yl)butan-2-yl)pyrimidine

FC1=CC(F)=C(C(CN2C=NC=N2)[C@@H](C)C3=CC=NC=N3)C=C1

20

4-(2-((2S,3R)-3-(2,4-difluorophenyl)-3-hydroxy-4-(1H-1,2,4-triazol-1-yl)butan-2-yl)-2,3-dihydrothiazol-4-yl)benzonitrile

FC1=C([C@](CN2N=CN=C2)(O)[C@H](C)C3NC(C4=CC=C(C#N)C=C4)=CS3)C=CC(F)=C1

 

Table 2: Molecular physicochemical descriptors analysis of 20 triazole analogues using Molinspiration online software tool.

Ligands

Log A a

TPSAb

Natomsc

MWd

nONe

nOH NH f

Nviolationsg

Nrotbh

Volumei

1

-0.6

41.58

5

69.07

3

1

0

0

60.72

2

0.90

30.72

9

125.17

3

0

0

3

128.06

3

2.92

30.72

13

181.28

3

0

0

7

195.27

4

-1.2

56.74

10

140.19

4

2

0

4

139.59

5

5.24

30.72

19

265.44

3

0

1

11

195.60

6

0.22

56.74

9

144.20

4

2

0

1

123.41

7

0.54

59.82

12

184.22

5

1

0

1

149.94

8

1.53

30.72

13

191.26

3

0

0

1

166.73

9

3.36

72.54

22

294.31

6

1

0

4

260.74

10

1.06

85.09

18

245.24

6

1

0

3

211.61

11

0.57

117.35

9

125.09

7

3

0

1

96.90

12

-0.1

81.66

22

306.28

7

1

0

5

248.96

13

2.0

95.38

25

343.39

9

0

0

6

310.42

14

-0.7

97.84

25

347.42

9

2

0

3

322.97

15

-2.4

143.11

20

301.20

10

3

0

5

235.61

16

-0.1

100.64

18

269.18

7

3

0

3

198.70

17

3.27

61.44

23

251.24

6

0

0

5

291.59

18

2.73

88.14

38

522.60

9

1

1

10

469.26

19

1.75

56.50

23

315.33

5

0

0

5

272.48

20

2.88

86.77

31

439.49

6

2

0

6

369.28

aOctanol-Water partition coefficient, b Polar surface area,  cNumber of non-hydrogen atoms,d Molecular weight, e Number of hydrogen bond acceptors [ O and N atoms], f Number of hydrogen bond donors [ OH and NH groups], g Number of Rule of 5 violations, h Number of rotatable bonds, i Molecular volume.

With regard to drug-likeness score, if the score is > 0 is active, -5.0 to -0.0 is moderate active and < -5.0 is inactive18.

 


It could be beneficial to know the physiochemical, drug-likeness and ADME (Absorption, Distribution, Metabolism and Excretion) properties of these 20 ligands before carry out docking studies. Lipinski’s rule of five was applied to know the above said properties and violation of the Lipinski’s rule of five is when logA>5, MW >500, number of N, O (hydrogen bond receptor) >10, number of OH and NH (hydrogen bond donor) >5 and number of the rotatable bond (rotb) >1517. In the present study, all the ligands (expect ligand 5 and 18) showed nil violation and complied well with the Lipinski’s rule of five as shown in Table 2.

All the ligands showed active to moderate active score towards all the six descriptions. Interestingly, none of them showed inactive score as shown in Table 3.


 

Table 3: Drug-likeness property analysis of 20 triazole analogues using Molinspiration tool.

Ligands

GPCR* ligand

Ion channel modulator

Kinase inhibitor

Nuclear receptor ligand

Protease inhibitor

Enzyme inhibitor

1

-3.72

-3.67

-3.68

-3.87

-3.82

-3.70

2

-2.06

-1.70

-2.53

-2.65

-2.51

-1.63

3

-0.52

-0.30

-0.85

-0.90

-0.85

-0.22

4

-0.56

-0.20

-0.85

-1.52

-0.84

-0.22

5

0.02

-0.35

-0.54

-0.34

-0.38

0.01

6

-3.60

-3.20

-3.04

-3.87

-3.68

-3.00

7

-1.74

-1.69

-1.55

-2.32

-1.58

-0.97

8

-1.35

-1.07

-1.58

-2.12

-1.73

-1.17

9

-0.67

-0.93

-0.50

-0.69

-0.84

-0.47

10

-0.40

-0.73

-0.98

-0.28

-0.65

-0.44

11

-1.90

-1.11

-1.88

-3.56

-2.48

-1.10

12

0.04

-0.01

-0.09

-0.01

-0.04

0.08

13

0.03

-0.12

-0.40

-0.18

-0.09

0.07

14

0.32

-0.30

0.07

-0.68

0.42

-10

15

0.59

0.67

0.12

-0.33

0.46

1.02

16

0.49

0.26

0.08

-0.81

-0.04

0.86

17

0.17

0.07

0.02

-0.29

-0.13

0.15

18

0.19

-0.11

0.01

-0.33

-0.04

0.10

19

0.25

0.01

0.06

-0.50

-0.02

0.04

20

0

-0.16

-0.23

-0.25

0.06

0.06

*GPCR- G Protein-coupled receptors

 

Table 4: ADME analysis of 20 triazole analogues using Swiss ADME online tool.

Ligands

Gl

BBB#

P-gp

CYP1A2*

CYP2C19*

CYP2C9*

CYP2D6*

CYP3A4*

Log Kp

1

High

No

No

No

No

No

No

No

-7.10

2

High

Yes

No

No

No

No

No

No

-6.63

3

High

Yes

No

Yes

No

No

No

No

-5.43

4

High

No

No

No

No

No

No

No

-7.77

5

High

Yes

No

Yes

Yes

No

Yes

No

-4.26

6

High

No

No

No

No

No

No

No

-6.83

7

High

No

No

No

No

No

No

No

-6.64

8

High

Yes

No

Yes

No

No

No

No

-6.07

9

High

Yes

No

Yes

Yes

No

No

No

-5.94

10

High

No

No

No

No

No

No

No

-6.48

11

High

No

No

No

No

No

No

No

-6.79

12

High

No

No

No

Yes

No

No

No

-8.62

13

High

No

No

No

Yes

No

No

No

-6.65

14

High

No

No

No

Yes

No

No

No

-8.15

15

Low

No

No

No

No

No

No

No

-10.1

16

High

No

No

No

No

No

No

No

-8.65

17

High

Yes

No

No

Yes

Yes

No

No

-6.17

18

High

Yes

Yes

Yes

Yes

Yes

Yes

Yes

-7.83

19

High

Yes

No

No

Yes

No

No

No

-6.46

20

High

No

Yes

No

Yes

Yes

Yes

Yes

-6.38

-Gastrointestinal, #-Blood-brain barrier permeant, -P-gp-P-glycoprotein substrate, *-CYP-Cytochrome P450 Inhibitors, -Skin Permeation (cm/s).

 


Absorption, Distribution, Metabolism and Excretion (ADME) is essential analysis tool and which is commonly accepted in the early stage of drug discovery, drug design and drug screening, owing to its unique characteristic nature19,20. Table 4 showed the ADME profile of the selected 20 triazole analogues; all the ligands are predicted to have high gastrointestinal (GI) absorption effect except ligand 15 [(2-(4-((2,4-dioxo-3,4-dihydropyrimidin-1(2H)-yl)methyl)-2H-1,2,3-triazol-2-yl)ethyl)phosphonic acid].

 

Human nucleophil elastase (HNE) is a key enzyme which plays major role in degenerative and inflammatory diseases, through proteolysis extracellular matrix (ECM) components21-24. The docking studies and binding site analyses in Table 5, shows that ligand 18 [2-(2,4-difluorophenyl)-1-(4-((1-(4-ethylbenzyl)-1H-1,2,3-triazol-4-yl)methyl)piperazin-1-yl)-3-(1H-1,2,4-triazol-1-yl)propan-2-ol] with the highest ACE [atomic contact energy (-330.87 kcal/mol)] while ligand 4 [4-(4H-1,2,4-triazol-4-yl)butan-1-amine] showed the least ACE (-79.11 kcal/mol) with that of HNE.

Table 5: The interaction energy analysis of 20 triazole analogues with Human nucleophil elastase (HNE) using PatchDock.

Ligands

-ACE

(-kcal/mol)

Interaction of amino acid residue

Bond distance (Å)

1

96.56

No Interaction

-

2

122

No Interaction

-

3

160.54

Ser 195

Ser 195

2.9

3.2

4

79.11

No Interaction

-

5

209.35

No Interaction

-

6

179.20

Cys 42 and Cys 42

His 57

2.3 and 2.5

3.4

7

211.03

Cys 58

3.5

8

197.29

No Interaction

-

9

241.84

His 57

2.0

10

228.60

No Interaction

-

11

127.67

Arg 129

3.0

12

225.13

No Interaction

-

13

248.58

His 57

3.1

14

273.56

Ser 195

1.9

15

182.89

Ser 214

2.5

16

116.31

No Interaction

-

17

267.03

Ser 195

2.7

18

330.87

Ser 195

3.5

19

182.51

Ser 195

2.6

20

276.56

No Interaction

-

 

The present finding was in good agreement with earlier report, where compound 3b [(N,N’,N,N’)-N,N’-(4,4’-(butane-1,4-diylbis(oxy))bis(4,1-phenylene)bis(methan-1-yl-1-ylidene)bis (3,5- diethyl-4H-1,2,4-triazol-4-amine)] exhibited potent elastase inhibitory activity25. In the present study, interaction with Ser 195 amino acid residue has been shown by five ligands such as ligand 3, 14, 17, 18 and 19.  Similarly, interaction with His 57 amino acid residue has been shown by three ligands such as ligand 6, 9, and 13. Previously Vergellia et al reported that Ser 195 amino acid present in the active site, (responsible for the nucleophilic attack) and similarly His 57 amino acid present in the catalytic site (in all the serine proteases)26.

 

Matrix metalloproteinase’s (MMPs) are family of zinc-dependent endopeptidases regulates both development and physiological events. There are now more than 20 members of the MMP family among them, MMP 2 (Gelatinase A) and MMP 9 (Gelatinase B) have been reported to be elevated in the pathological conditions like as aging, cancer, inflammation and wound healing21.  The docking studies and binding site analyses in Table 6, shows that ligand 18 [2-(2,4-difluorophenyl)-1-(4-((1-(4-ethylbenzyl)-1H-1,2,3-triazol-4-yl)methyl)piperazin-1-yl)-3-(1H-1,2,4-triazol-1-yl)propan-2-ol] with the highest ACE [atomic contact energy (-490.13 kcal/mol)] while ligand 1 [4H-1,2,4-triazole] showed the least ACE (-49.60 kcal/mol) with that of MMP 2. Similarly in the case of MMP 9, Table 7 shows that ligand 18 [2-(2,4-difluorophenyl)-1-(4-((1-(4-ethylbenzyl)-1H-1,2,3-triazol-4-yl)methyl)piperazin-1-yl)-3-(1H-1,2,4-triazol-1-yl)propan-2-ol] with the highest ACE [atomic contact energy (-489.56 kcal/mol)] while ligand 1 [4H-1,2,4-triazole] showed the least ACE (-112.43 kcal/mol) with that of MMP 9.

 

Table 6: The interaction energy analysis of 20 triazole analogues with MMP 2 using PatchDock

Ligands

-ACE

(-kcal/mol)

Interaction of amino acid residue

Bond distance (Å)

1

49.60

Leu 185

3.2, 3.6

2

181.45

No Interaction

-

3

207.11

No Interaction

-

4

190.27

No Interaction

-

5

317.72

No Interaction

-

6

188.75

Pro 221

Ala 220

2.9

3.4

7

227.17

Thr 229 and Leu 218

Ala 217

2.9 and 2.9

3.0

8

274.32

No Interaction

-

9

378.02

No Interaction

-

10

273.05

No Interaction

-

11

195.66

Ile 222 and Ala 217

3.5 and 3.4

12

303.82

Thr 229 and Arg 233

2.2 and 3.6

13

335.88

Ile 222

3.0

14

391.27

Glu 202

2.4, 2.6

15

107.73

Ala 167 and Ala 165

3.1 and 3.5

16

176.82

Thr 227

3.3

17

357.87

No Interaction

-

18

490.13

Leu 164 and Ala 165

2.3 and 2.9

19

289.68

Thr 229

2.9

20

377.11

No Interaction

-

 

Table 7:  The interaction energy analysis of 20 triazole analogues with MMP 9 using PatchDock.

Ligands

-ACE

(-kcal/mol)

Interaction of

amino acid residue

Bond distance (Å)

1

112.43

No Interaction

-

2

166.06

No Interaction

-

3

135.59

Pro 240 and Ala 242

2.9 and 2.9

4

200.30

No Interaction

-

5

326.86

No Interaction

-

6

235.67

Met 247 and Ala 242

2.9 and 3.3

7

269.80

Ala 242 and Leu 243

3.0 and 3.4

8

315.72

Pro 240

2.7

9

375.31

No Interaction

-

10

303.74

No Interaction

-

11

195.12

Arg 249 and Leu 222

3.3 and 3.4

12

151.73

His 230 and Gln 227

3.2 and 3.5

13

284.94

No Interaction

-

14

398.42

Met 247 and Leu 243

Pro 240 and His 266

Tyr

3.0 and 2.7, 3.2

3.4 and 3.2

2.9

15

243.79

Arg 249

3.1, 3.4

16

180.57

No Interaction

-

17

193.44

No Interaction

-

18

489.56

Leu 188

3.5

19

310.21

Pro 254 and Pro 240

3.3 and 2.8

20

416.42

Leu 188

3.1

 

Triazole derivatives containing thiirane compound reported to exhibit strong gelatinase inhibitory activity28. In the present study, interaction with Leu 188 amino acid residue has been shown by two ligands such as ligand 18 [2-(2,4-difluorophenyl)-1-(4-((1-(4-ethylbenzyl)-1H-1,2,3-triazol-4-yl)methyl)piperazin-1-yl)-3-(1H-1,2,4-triazol-1-yl)propan-2-ol] and ligand 20 [4-(2-((2S,3R)-3-(2,4-difluorophenyl)-3-hydroxy-4-(1H-1,2,4-triazol-1-yl)butan-2-yl)-2,3-dihydrothiazol-4-yl)benzonitrile]. The present finding was in good agreement with earlier report21.

 

Tyrosinase is the main regulatory enzyme in melanogenesis pathway that too mainly in the primary two stages such as (i) tyrosine hydroxylation to 3, 4-dihydroxyphenylalanine (DOPA) and (ii) the oxidation of DOPA to dopoquinone23.  The docking studies and binding site analyses in Table 8, shows that ligand 10 [2-(4-acetylphenyl)-5-methyl-2H-1, 2, 3-triazole-4-carboxylic acid] with the highest ACE [atomic contact energy (-152.56 kcal/mol)]. While four ligands such as ligand 12 [2-(2,4-difluorophenyl)-1,3-di(4H-1,2,4-triazol-4-yl)propan-2-ol], ligand 16 [(2R,3S,4R,5R)-2-(hydroxymethyl)-5-(5-(trifluoromethyl)-1H-1,2,3-triazol-1-yl)tetrahydrofuran-3,4-diol], ligand 19 [4-((2R)-3-(2,4-difluorophenyl)-4-(1H-1,2,4-triazol-1-yl)butan-2-yl)pyrimidine] and ligand 20 [[4-(2-((2S,3R)-3-(2,4-difluorophenyl)-3-hydroxy-4-(1H-1,2,4-triazol-1-yl)butan-2-yl)-2,3-dihydrothiazol-4-yl)benzonitrile] showed poor binding phenomenon as reported by Akdogan and co-workers29.

 

Table 8: The interaction energy analysis of 20 triazole analogues with tyrosinase inhibitors using PatchDock.

Ligands

-ACE

(-kcal/mol)

Interaction of amino acid residue

Bond distance (Å)

1

16.16

Ser 364

3.3, 3.4

2

51.89

No Interaction

-

3

47.54

No Interaction

-

4

0.96

Gln 74

2.8

5

35.19

No Interaction

-

6

112.16

Gln 307

2.9

7

87.08

No Interaction         

-

8

117.45

Lys 379

3.2, 3.5

9

95.90

Thr 308

2.7, 3.3

10

152.56

No Interaction         

-

11

82.05

Asn 243

Gly 86

2.6

3.4

12

+ 27.80

Gln 307

2.5, 3.2

13

72.40

Asp 312

3.1

14

54.32

Tyr 82

3.2

15

39.09

Lys 379

Thr 308

Gln 307

1.9, 2.8

2.9

3.5

16

+ 36.02

Asp 312

2.8

17

60.77

Lys 379

2.6

18

115.92

No Interaction

-

19

+ 41.46

Asp 312

3.3

20

+ 18.49

Lys 372

3.2

 

2,2’-(4,4’-(butane-1,4-,diyl) bis (3-methy1-5-oxo-1H-1,2,4-triazole-4,1(4H,5H)-diyl))bis(N’-(thiophen-2-ylmethylene)acetohydrazide and 2,2’-(4,4’-(hexane-1,6-,diyl) bis (3-methy1-5-oxo-1H-1,2,4-triazole-4,1(4H,5H)-diyl))bis(N’-(thiophen-2-ylmethylene)acetohydrazide have been reported to inhibit tyrosinase activity30. Similarly, (Z)-3-(2-Fluorobenzyl)-4-[(4-fluorobenzylidene) amino]-1H-1, 2, 4-triazole-5(4H)-thione and (Z)-3-(2-Fluorobenzyl)-4-[(4-hydroxybenzylidene) amino]-1H-1, 2, 4-triazole-5(4H)-thione have been reported to exhibit strong tyrosinase inhibitory activity31. In the present study, interaction with Lys 379 amino acid residue has been shown by three ligands such as ligand 8 [5-methyl-4-phenyl-4H-1,2,4-triazole-3-thiol], ligand 15 [(2-(4-((2,4-dioxo-3,4-dihydropyrimidin-1(2H)-yl)methyl)-2H-1,2,3-triazol-2-yl)ethyl)phosphonic acid] and ligand 17 [1-(1-(2,4-dichlorophenyl)-2-(1H-1,2,3-triazol-1-yl)ethyl)-4-ethyl-5-methyl-1H-1,2,3-triazole]. The present finding was in good agreement with earlier report23.

 

CONCLUSION:

In the present study ligand 18 [2-(2,4-difluorophenyl)-1-(4-((1-(4-ethylbenzyl)-1H-1,2,3-triazol-4-yl)methyl) piperazin-1-yl)-3-(1H-1,2,4-triazol-1-yl) propan-2-ol] had exhibited the maximum binding energy for three target enzymes such as HNE, MMP 2 and MMP 9, whereas in the case of tyrosinase ligand 10 [2-(4-acetylphenyl)-5-methyl-2H-1, 2, 3-triazole-4-carboxylic acid] had shown the maximum binding energy. In addition to this all the ligands (expect ligand 5 and 18) exhibited zero violation and complied very well with the Lipinski’s rule of five. And moreover, ADME analysis showed that all the ligands (expect ligand 15) predicated to have high gastrointestinal (GI) absorption effect. Hence, the present study has paved a new insight in understanding 20 triazole analogues as potential inhibitors against HNE, MMP 2, MMP 9 and tyrosinase.

 

ACKNOWLEDGEMENTS:

The author VV thankful to Department of Chemistry, Anna University, Chennai-25 and N.R would like to express his sincere gratitude to Department of Biochemistry, St. Peter’s Institute Higher Education and Research (Deemed to be University), Avadi, Chennai-54, Tamil Nadu, India for providing research facilities.

 

CONFLICT OF INTEREST:

The authors declare no conflict of interest.

 

REFERENCES:

1.      Asif M. A brief review on anti-tubercular activity of pharmacological activity of some triazole analogues. Global J Res and Rev. 2014; 1(3):051-058.

2.      Otuokere IE, Amaku FJ, Alisa CO. In silico geometry optimization, excited – state properties of (2E)-N-Hydroxy-3-[3-(Phenylsulfamoyl) Phenyl] prop-2-enamide (Belinostat) and its molecular docking studies with Ebola Virus glycoprotein. J Pharm Res 2015; 5(3):131-137.

3.      Sudhakar P, Pushkalai S, Poorana, Sabarinath C, Priyadharshini S, Haripriya S. Molecular docking and synthesis of 1, 2, 4 -triazin analogue of diclofenac as potential ligand for parkinson's. Res J Pharm  Pharmcodynamics 2018; 10(1):8-12.

4.      Panigrahi N, Ganguly S, Jagadeesh P. Synthesis, antimicrobial evaluation and molecular docking studies of novel Oxazolidinone-thiophene chalcone hybrid derivatives. A J Res and Tech 2018; 11(12):5611-5622.

5.      Chellamalla R, Makula A. Molecular docking studies and ADMET predictions of pyrimidine coumarin scaffolds as potential IDO inhibitors. A J Res Chem 2017; 10(3):331-340.

6.      Hasan T, Ghalib RM, Mehdi SH, Singh PK, Baqri SSR. Normal mode analysis, electronic parameters and molecular docking study of 3, 5, 4’-Trihydroxy-6, 7-dimethoxy-flavone (Eupalitin) using first principle. A J Res Chem 2017; 10(6):789-797.

7.      Ramjith US, Shahim M. Molecular docking study of novel imidazo-1, 3, 4 thiadiazole derivatives. Res J Pharm and Tech 2013; 6(6):688-694.

8.      Pavlo V, Zadorozhnii VV, Kiselev, AE, Titova, AV, Kharchenko, Ihor O. Pokotylo OV, Okhtina. Molecular docking studies of N-5-aryl-1, 3, 4-oxadiazolo-2, 2-dichloroacet-amidines as inhibitors of enoyl-ACP reductase Mycobacterium tuberculosis. Res J Pharm and Tech 2017; 10(4):1091-1097.

9.      Pavlo VZadorozhnii, Vadym VK, Nataliia OT, Aleksandr VK, Ihor OP, Oxana VO, Oxana VK. In Silico prediction and molecular docking studies of N-Amidoalkylated derivatives of 1, 3, 4-oxadiazole as COX-1 and COX-2 potential inhibitors. Res J Pharm and Tech 2017; 10(11):3957-3963.

10.   Godhani DR, Jogel AA, Sanghani AM, Mehta JP. Synthesis and biological screening of 1, 2, 4-triazole derivatives. Ind J Chem. 2015; 54B:556-564.

11.   Kumar H, Javed SA, Khan SA, Amir M. 1, 3, 4-Oxadiazole/thiazole and 1,2,4-triazole derivatives of biphenyl-4-yloxy acetic acid: Synthesis and preliminary evaluation of biological properties. Euro J Med Chem. 2008; 43:2688-2698.

12.   Wuest F, Tang X, Kniess T, Pietzsch J, Suresh M. Synthesis and cyclooxygenase inhibition of various (aryl-1,2,3-triazole-1-yl)-methanesulfonylphenyl derivatives. Bio Med Chem. 2009; 17:1146-1151.

13.   Vijayakumar V, Radhakrishnan N, Rameshkumar C. Molecular docking analysis of imidazole derivatives and polybenzimidazole analogues as inhibitors of superoxide dismutase (SOD) and Xanthine oxidase (XO). IEEE International Conference on Smart Technologies and Management for Computing, Communication, Controls, Energy and Materials (ICSTM) 2017; 513-516.

14.   Navjot Kaur, Monika, Kulwinder Singh. 3D-QSAR and Molecular Docking Studies of N-(2-Aminophenyl)-Benzamide Derivatives as Inhibitors of HDAC2. Res J Pharm and Tech 2014; 7(7):760-770.

15.   Wang J, Urban L. The impact of early ADME profiling on drug discovery and development strategy. Drug Discovery World. 2004; 5:73-96.

16.   Hailu M, Narayanaswamy R, Abrehame S, Manoj VR. In silico analysis of Glycogen Synthase Kinase 3-Beta (GSK-3β) as aflatoxin binder. Biomed Pharmacol J 2019; 12(3):1449-1456.

17.   Lipinski CA, Lambardo F, Dominy BW, Feeney PJ. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev. 2001; 46:3-26.

18.   Singh S, Gupta AK, Verma A. Molecular properties and bioactivity score of the Aloe vera antioxidant compounds- In order to lead finding. Res J Pharm Biol Chem Sci. 2013; 4(2):876-891.

19.   Vijayakumar V, Radhakrishnan N, Rameshkumar C. Predicting the biodegradability  nature of imidazole and its derivatives by modelling two histidine degradation enzyme (Urocanase and Formiminoglutamase) activities. Asian J Pharm Clin Res. 2017; 10 (11):383-386.

20.   Wang J, Urban L. The impact of early ADME profiling on drug discovery and development strategy. Drug Discovery World. 2003; 5:73-96.

21.   Radhakrishnan N, Lam KW, Intan IS. In silico analysis of selected honey constituents as human neutrophil elastase (HNE) and matrix metalloproteinases (MMP 2 and 9) inhibitors. Int J Food Prop. 2015; 18:2155-2164.

22.   Radhakrishnan N, Lam KW, Faridah A, Intan IS. Molecular docking analysis of curcumin analogues elastase (HNE) inhibitors. Bangladesh J. Pharmacol. 2014; 9:77-82.

23.   Radhakrishnan N, Lam Kok Wai, Norhaizan Mohd Esa. Molecular docking analysis of phytic acid and 4-hydroxyisoleucins as cyclooxygenase-2, microsomal prostaglandin E synthase-2, tyrosinase, human neutrophil elastase, matrix metalloproteinase-2 and 9, xanthine oxidase, squalene synthase, nitric oxide synthase, human aldose reductase and lipoxygenase inhibitors. Pharmacogn Mag. 2017; 13 (51):S512-S518.

24.   Vijayakumar V, Prabakaran A, Radhakrishnan N, Muthu S, Rameshkumar C, Paulraj EI. Synthesis, characterization, spectroscopic studies, DFT and molecular docking analysis of N4, N4′-dibutyl-3, 3′-diaminobenzidine. J Mol Struct. 2019; 1179:325-335.

25.   Gumrukcuoglu N, Sokmen BB, Ugras S, Ugras HI, Yanardag R. Synthesis, antibacterial, antielastase, antiurease and antioxidane activities of new 1,4-butylene bridged bis-1,2,4-triazole derivatives. J Enzyme Inhib Med Chem. 2013; 28(1):89-94.

26.   Vergelli C, Schepetkin IA, Crocetti L, Lacovone A, Giovannoni MP, Guerrini G, Khlebnikov AI, Ciattini S, Ciciani G, Quinn MT. Isoxazol-5(2H)-one: a new scafford for potent heman neutrophil elastase (HNE) inhibitors. J Enzyme Inhib Med Chem. 2017; 32(1):821-831.

27.   Romanchikova N, Trapencieris P, Zemitis J, Turks M. A new matrix metalloproteinase-2 inhibitor triazolylmethyl aziridine reduces melanoma cell invasion, angiogenesis and tergets ERK1/2 phosphorylation. J Enzyme Inhib Med Chem. 2014; 29(6):765-772.

28.   Testero SA, Llarrull LI, Fisher JF, Chang M, Mobashey S. Exploring the functional space of thiiranes as gelatinase inhibitors using click chemistry. ARKIVOC 2011; 7:221-236.

29.   Akdogan ED, Erman B, Yelekci K. In silico design of novel and highly selective lysine‑specific histone demethylase inhibitors. Turk J Chem. 2011; 35:523‑42.

30.   Unver Y, Celik F, Dugdu E, Barut B, Sancak K. Synthesis of new bis-triazole compounds including Schiff bases and enzyme inhibition and antioxidant activities. J  Med Chem. 2016; 1(1): 005-010.

31.   Rafiq M, Saleem M, Hanif M, Kang SK, Seo SY, Lee KH. Synthesis, structural elucidation and bioevaluation of 4-amino-1, 2, 4-triazole-3-thione’s Schiff base derivatives. Arch Pharm Res. 2015. 

 

 

 

Received on 13.09.2019            Modified on 28.10.2019

Accepted on 14.12.2019           © RJPT All right reserved

Research J. Pharm. and Tech 2020; 13(6): 2777-2783.

DOI: 10.5958/0974-360X.2020.00493.X